METHOD FOR PREDICTING AN AMOUNT OF DIMENSIONAL ACCURACY 
DEFECT AT THE TIME OF PRESS-FORMING METAL SHEET 

Background of the Invention 

1. Field of the Invention 

The present invention relates to a method for predicting 
an amount of dimensional accuracy defect which is caused when 
a metal sheet such as a thin steel sheet, an aluminum sheet or 
the like which is mainly applied to an automobile body is formed 
by a press, and more particularly to a method which can 
preliminarily simply and accurately predict an amount of 
dimensional accuracy defect (mainly, a wall warp amount, an 
angular change amount or the like) of a formed product caused 
by an elastic recovery after removing the product from a mold 
in a press forming. 

2 . Description of the Prior Art 

Many parts of an automobile body are generally constituted 
by parts which are produced by press-forming thin steel sheets. 
However, in forming these parts by a press-forming, due to an 
elastic recovery behavior after the removal from mold (the 
removal of a formed product from a mold), the shape (dimension) 
of the formed product is changed from the designed value and there 
exists a case that drawbacks take place at the time of assembling 
or joining formed products. The joining is mainly a joining by 
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a spot welding. These drawbacks are generally called 
dimensional accuracy defects. As such dimensional accuracy 
defects, the wall warp and the angular change or the like are 
known. 

Recently, from the viewpoint of making the automobile body 
light-weighted and increasing the stability of the automobile, 
cases that thin steel sheets having a high strength or aluminum 
sheets being light-weighted but having a lower Young's modulus 
compared to the steel sheets are used as the material for the 
automobile body have been increasing and hence, the above- 
mentioned dimensional accuracy defects have been taken as 
crucial problems. 

To solve such problems, as means for enhancing the 
dimensional accuracy, (A) the development of new techniques and 
(B) the adjustment by predictive techniques and the like have 
been mainly carried out. As a technique belonging to (A), 
inventors of the present invention have proposed a method which 
reduces a wall warp amount by applying a tension to wall portion 
by increasing the blank holding force at the final state of the 
press operation. As a technique belonging to (B) , there has been 
known a method which predicts an elastic recovery amount after 
the removal of the sheet from a mold and determines the dimension 
of a mold with a value which is obtained by subtracting the 
predicted value of the elastic recovery amount from the target 
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product dimension and then adjusts the dimension to the 
prescribed dimension after removing the sheet from the mold. 

With respect to the above-mentioned two means, whichever 
means is adopted, it is extremely important to preliminarily 
accurately predict an amount of dimensional accuracy defect 
after forming based on information on a metal sheet to be formed, 
forming conditions and the like. 

Heretofore, in predicting an amount of dimensional 
accuracy defect after press-forming a metal sheet, there have 
been known (1) a method which performs the prediction based on 
the experience of skilled workers, the stored past results and 
the like, (2) a method which performs the prediction based on 
a numerical simulations and the like. However, it has been 
pointed out that either one of these methods suffers from 
following problems. 

With respect to the method (1), it is difficult for 
engineers having a little experience or a little storage of 
techniques to perform an accurate prediction and hence, the 
method (1) is not suitable for a practical operation. On the 
other hand, with respect to the method (2), although various 
techniques which are proposed in "Sosei to Kako" (for example, 
vol. 36, No. 410 (1995) p203 - 210, vol.37. No. 410(1996) pl352 
- 133 6 and the like) have been known, the expertise such as 
numerical simulations, mathematics or the like is necessary and 

-3- 



there also exists a problem that a simulation facility such as 
a computer becomes large-sized. 

Summary of the Invention 
The present invention has been made under the above- 
mentioned background and it is an object of the present invention 
to provide a method which enables engineers who do not have the 
experience and the storage of technique and also do not have the 
expertise such as the numerical simulation, mathematics and the 
like can preliminarily simply and accurately predict an amount 
of dimensional accuracy defect at the time of press-forming a 
metal sheet. 

The gist of a prediction method of the present invention 
capable of achieving the above-mentioned object lies in that, 
in predicting an amount of dimensional accuracy defect at the 
time of press-forming a metal sheet, as a stress-strain 
relationship, an elastic-perf ectly plastic solid model having 
a fixed stress value is adopted after yielding and a value which 
is equal to or less than the tensile strength and exceeds the 
yield strength is called an apparent yield strength. 

As the specific constitution of the above-mentioned 
method, a constitution which uses an internal division value 
obtained by a following equation (1) based on the yield strength 
(YS) and the tensile strength (TS) which are actually measured 
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as the above-mentioned apparent yield strength ( fJ p ' ) is 
considered. 

Cp' = k • YS + (1 - k)TS (1) 
where k is a coefficient of the internal division value 

Further, the coefficient k in the above-mentioned 
equation (1) can be expressed as a function [k = f (TS/t) ] of the 
ratio (TS/t) between the tensile strength TS and the plate 
thickness t and takes a value in the range of 0 < k < 1 . 

In the above-mentioned method of the present invention, 
in predicting a wall warp amount at the time of press-forming, 
the wall warp amount is predicted based on p which is obtained 
by following equations (2) - (4). 

<3= (3 O-p' /Et) • < 1 - D • [ ( CTt/TS) - 0.3]^} 

- C • ( rd - 5 ) ( 2 ) 

crp' = k • YS + (1 - k)TS (3) 

k = A • (TS/t) + B (4) 
where p = wall warp amount (curvature; l/mm),crp': apparent 
yield strength (MPa), E: Young's modulus (MPa), t: sheet 
thickness (mm) , TS : tensile strength (actually measured value; 
MPa), YS: yield strength (actually measured value; MPa),0'T: 
tension acting on wall portion (MPa), rd: radius of die shoulder 
of press forming tool (mm), k: coefficient of internal division 
value. A: negative constant, B, C, D: positive constant. 

Further, in the above-mentioned method of the present 
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invention, in predicting an angular change amount at the time 
of press-forming, the angular change amount is predicted based 
on A 0 which is obtained by following equations (5) - (8). 
A 0= • (rp + t/2) • A p (5) 
A p= (-3 (7p'/Et) • [1 + exp(-G • rp) ] (6) 
crp' = k • YS + ( 1 - k)TS (7) 
k = A • (TS/t) + B (8) 
where A^: angular change amount (degree), 6: bending angle 
(degree), rp: radius of shoulder of bending tool (mm), t: sheet 
thickness (mm), A p : curvature change amount (1/mm), crp': 
apparent yield strength (MPa), E: Young's modulus (MPa), YS: 
yield strength (actually measured value; MPa),TS: tensile 
strength (actually measured value; MPa), k: coefficient of 
internal division value. A: negative constant, B, G: positive 
constant . 

Brief Description of the Drawings 
Fig. 1 is an explanatory view showing an example of an outer 
appearance of a "hat channel" specimen. 

Fig. 2 is a schematic explanatory view showing a main 
forming method of a "hat channel" specimen. 

Fig. 3 is a view for explaining the dimensional accuracy 
defects which are subjects of the present invention. 

Fig. 4 is graphs showing the relationship between major 
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influence factors which affect a wall warp and a wall warp amount 
P . 

Fig. 5 is a stress - strain diagram used in a conventional 
prediction method. 

Fig. 6 is a stress - strain diagram when an elastic- 
perfectly plastic solid material (work hardening not considered) 
is assumed. 

Fig. 7 is a view for explaining the deformation history 
of a wall portion of material subjected to a press-forming. 

Fig. 8 is a view showing the distribution of the strain 
e and the stress a acting in the plate thickness direction of 
the material (at the marked point in Fig. 7) and the distribution 
of an elastic region and a plastic region in the inside of the 
material at a stage of Fig. 7(a). 

Fig. 9 is a view showing the distribution of the strain 
£ and the stress a acting in the plate thickness direction of 
the material (at the marked point in Fig. 7) and the distribution 
of an elastic region and a plastic region in the inside of the 
material at a stage of Fig. 7(b). 

Fig. 10 is a view showing the distribution of the strain 
£ and the stress O acting in the plate thickness direction of 
the material (at the marked point in Fig. 7) and the distribution 
of an elastic region and a plastic region in the inside of the 
material at a stage of Fig- 7(c). 
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Fig. 11 is a view showing the distribution of the strain 
£ and the stress a acting in the plate thickness direction of 
the material (at the marked point in Fig. 7) and the distribution 
of an elastic region and a plastic region in the inside of the 
material after removing a product from a mold. 

Fig. 12 is a stress - strain diagram used in the present 
invention. 

Fig. 13 is a graph showing the relationship between major 
influence factors which affect an angular change and an angular 
change amount A Q . 

Fig. 14 is a graph showing in comparison an actually 
measured value (pmes) and a predicted value (pcal) based on 
the above-mentioned equations (2) - (4) with respect to wall warp 
amounts (curvatures) p of various kinds of steel sheets shown 
in Table 1 . 

Fig. 15 is a graph showing in comparison an actually 
measured value (pmes) and a predicted value (pcal) based on 
the above-mentioned equations (2) - (4) with respect to wall warp 
amounts (curvatures) p of various kinds of aluminum sheets 
shown in Table 2. 

Fig. 16 is a graph showing in comparison an actually 
measured value (pmes) and a predicted value (peal) based on 
the above-mentioned equations (2) - (4) with respect to wall warp 
amounts (curvatures) p when the tension (Jt is relatively large . 
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Fig. 17 is a schematic view showing the state of an L- 
bend forming experiment. 

Fig. 18 is a graph showing in comparison an actually 
measured value (A^mes) and a predicted value (A^cal) based 
on the above-mentioned equations ( 5 ) - ( 8 ) with respect to angular 
change amounts A 0 of various kinds of steel sheets shown in Table 
3. 

Fig. 19 is a schematic view showing the state of a u- 
bend forming experiment. 

Fig. 2 0 is a graph showing in comparison an actually 
measured value (A^mes) and a predicted value (A^cal) based 
on the above-mentioned equations ( 5 ) - ( 8 ) with respect to angular 
change amounts A 0 at the time of U-bend forming. 

Fig. 21 is views showing the influence which the tensile 
strength and the thickness of sheet give to the work hardening. 

Description of the Preferred Embodiment 
Inventors of the present invention have carried out 
extensive studies from various aspects to solve the above- 
mentioned drawbacks. Eventually the inventors have made a 
finding that with the adoption of the above-mentioned 
constitution, the above-mentioned object of the present 
invention can be achieved. Following the procedure which has 
led to the completion of the present invention, the manner of 
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operation and the advantageous effects of the present invention 
are explained hereinafter in conjunction with attached drawings . 
In the explanation described hereinafter, for the explanation 
purpose, a case in which a "hat channel" specimen which is 
popularly used for parts of an automobile is formed as a material 
to be subjected to the prediction of an amount of dimensional 
accuracy defect is picked up and explained. However, it is 
needless to say that the material to be subjected to the 
prediction of an amount of dimensional accuracy defect in the 
present invention is not limited to such a "hat channel" specimen. 

Fig. 1 is an explanatory view showing an example of an 
outer configuration of the above-mentioned "hat channel" 
specimen. As a major method for forming such a "hat channel" 
specimen, as shown in Fig. 2, a drawing forming method (Fig. 2(a)) 
and a bending forming method (Fig. 2(b)) are known. In these 
forming methods, the dimensional accuracy defect which 
constitutes problem in particular is comprised of a "wall warp" 
phenomenon which occurs mainly at the time of performing a deep 
drawing and an "angular change" phenomenon which occurs at the 
time of performing a bending. 

Assuming that a design (target) shape (a cross-sectional 
shape perpendicular to an axis) of the "hat channel" specimen 
is a shape shown in Fig. 3(a), the "wall warp" phenomenon is a 
phenomenon in which a wall portion between R stops (R being the 
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radius of a die shoulder of a press forming tool) is warped as 
indicated by a broken line portion of Fig. 3(b). Further, 
assuming a target angle of a bending portion as 0 (Fig 3(a)), 
the "angular change" phenomenon is a phenomenon in which the 
bending portion is formed with an angle 0^ which is larger than 
the target angled as indicated by a broken line portion of Fig. 
3(c). The wall warp amount p is expressed by a curvature ( 1 /mm) 
of the above-mentioned R, while the angular change amount A 
6 is expressed by the difference between target angle 0 and the 
angle 0^ made after forming the hat channel" specimen. Then, 
the specific method for predicting these wall warp amount p and 
angular change amount A ^ is explained hereinafter in detail. 

At the time of performing the deep drawing, a metal sheet 
subjected to such a forming receives the bending/unbending 
deformation when the metal sheet passes a shoulder portion of 
a press forming tool (die shoulder). Accordingly, at the time 
of forming, in the wall portion, the difference between stresses 
having opposite signs (the difference between a tensile stress 
and a compression stress) is generated on the front and back in 
the sheet thickness direction (see Fig. 9 and Fig. 10 which will 
be explained later). Although a bending moment is generated in 
the sheet thickness direction due to this difference between 
stresses having opposite signs, when an external force is removed 
at the time of removing the metal sheet from the mold, the 
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deformation is recovered to some extent due to the elasticity 
which the material holds (elastic recovery). It is known that 
this elastic recovery behavior is the major cause of the wall 
warp. The elastic recovery behavior is generated so as to make 
the bending moment zero and the strain and the stress remain in 
the sheet thickness direction as shown in Fig. 11 (b), (c) which 
will be explained later. 

Further, it has been known that following tendencies (1) 
to (5) exist with respect to the wall warp amount p. Here, the 
relationship between the major influence factors which affect 
the wall warp and the wall warp amount p is shown in Fig. 4. 

( 1 ) Corresponding to the increase of the strength (tensile 
strength TS) of the metal material, the wall warp amount p is 
also increased (Fig. 4(a)). 

(2) Corresponding to the decrease of the sheet thickness 
t of the metal sheet, the wall warp amount p is increased (Fig. 
4(b)). 

(3) Corresponding to the decrease of the Young's modulus 
E of the metal material, the wall warp amount p is increased 
(Fig. 4(c)). 

(4) Corresponding to the increase of the radius rd of the 
die shoulder of the press forming tool, the wall warp amount p 
is increased (Fig. 4(d)). However, the region where the radius 
rd of the die shoulder of the press forming tool is extremely 
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small is excluded. 

(5) Corresponding to the increase of the tension (7i (blank 
holding force) working on a wall portion, the wall warp amount 
p is increased (a region indicated by a broken line in Fig. 4(e)). 
However, the region where the tension is extremely small is 
excluded. 

The inventors of the present invention have made extensive 
studies from various aspects aiming at the realization of a 
prediction equation which can accurately (quantitatively) and 
easily predict these tendencies based on the above-mentioned 
finding. To this end, first of all, the inventors of the present 
invention have reviewed reasons why the complicate calculations 
such as numerical value simulations have been carried out in the 
prior art to predict the wall warp amount. As a result, we 
arrived at a conclusion that one of the main reasons is that, 
in the prior art, as the stress-strain relationship at the time 
of applying the deformation (strain), the relationship which is 
close to the reality as shown in Fig. 5, that is, the work 
hardening in the plastic deformation range where the stress is 
increased corresponding to the increase of the plastic strain 
is taken into consideration as the direct subject of the numerical 
value simulation. 

Accordingly, to simplify the prediction equation, the 
inventors of the present invention have assumed the stress- 
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strain relationship that the metal material does not generate 
the work hardening although the relationship is different from 
the actual material behavior and then have tried to formulate 
the more simplified prediction equation under such a stress- 
strain relationship. That is, as shown in Fig. 6, the inventors 
have assumed a plastic material whose stress keeps a fixed value 
after yielding (a yielding stress (7p) {such a plastic material 
being generally called an "elastic-perf ectly plastic solid 
material") and then have tried to formulate the prediction 
equation based on such a plastic body. 

As a result, the inventors of the present invention have 
found that an amount of curvature change due to the elastic 
recovery after bending/unbending deformation (the wall warp 
amount p when the tension is extremely small) can be expressed 
by a following simple equation (I). Further, the inventors of 
the present invention also have found that when a die having a 
die shoulder diameter of 3 - 20 mm which is popularly used in 
a usual press forming is used, the equation (I) can be expressed 
by a further simplified equation such as a following equation 
(II). 

p = (3crp/Et) • { 1 - 7/3 • [(2crp/Et) • rd]^> (i) 
p = (3crp/Et) (II) 
where p = wall warp amount (curvature; 1/mm) , cTp: yield strength 
(MPa), E: Young's modulus (MPa), t: sheet thickness (mm) , rd: 
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radius of die shoulder of press forming tool (mm) . 

The above-mentioned equations (I) (II) can be formulated 
by analyzing the change of bending moment in the deformation steps 
at the time of press-forming. This procedure is explained in 
conjunction with drawings. Fig. 7 is a view for explaining the 
deformation history of a wall portion of a metal sheet subjected 
to the press-forming (particularly deep drawing). In the 
drawing, the metal sheet is sequentially formed by a press in 
the order of Fig. 7(a) -> Fig. 7(b) -> Fig. 7(c). A portion of 
the metal sheet where the wall warp occurs is indicated by a mark 
X as at the marked point in Fig. 7. Fig. 7(a) shows the stage 
before press-forming. Fig. 7(b) shows the bending stage, and Fig. 
7(c) shows the unbending stage. 

On the other hand. Fig. 8 to Fig. 10 show the distribution 
of the strain £ and the stress C acting in the sheet thickness 
direction of the metal sheet (the above-mentioned marked point 
X ) and the distribution of an elastic region and a plastic region 
in the inside of the metal sheet at respective stages when the 
metal sheet is formed by the press in the order of Fig. 7(a) -> 
Fig. 7(b) -> Fig. 7(c). Fig. 8 corresponds to the stage shown 
in Fig. 7(a), Fig. 9 corresponds to the stage shown in Fig. 7(b), 
and Fig. 10 corresponds to the stage shown in Fig. 7(c). The 
distribution of the strain £ and the stress cr acting in the 
plate thickness direction of the metal sheet (the above- 

-15- 



mentioned marked point X ) and the distribution of an elastic 
region and a plastic region in the inside of the metal sheet after 
removing the product from the mold are shown in Fig 11. 
Respective views (a) in Figs. 8 to 11 indicate the distribution 
of the elastic region and the plastic region in the inside of 
the metal sheet. Respective views (b) in Figs. 8 to 11 indicate 
the distribution of the strain £ . Respective views (c) in Figs. 
8 to 11 indicate the distribution of the stress a . 

First of all, in the stage shown in Fig. 7(a), the marked 
point X is formed of the elastic region over the entire region 
in the sheet thickness direction (Fig. 8(a)) and neither the 
distribution of the strains nor the distribution of the stress 
a is generated (Fig. 8(b), (c)). As a result, no bending moment 
is generated in this stage. Accordingly, the bending moment 
amount M® acting on the central portion of the sheet thickness 
t can be expressed by a following equation (III), wherein the 
position V in the sheet thickness direction is set as a variable. 

M^^^.o-rjcLri^^O (HI) 

2 

Then, in the stage (bending stage) shown in Fig. 7(b), 
in the marked point X , both end surface sides in the sheet 
thickness direction of the metal sheet become the plastic region 
and the central portion becomes the elastic region (Fig. 9(a)). 
Here, the distribution of the strain £ is made such that the 
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strains having the maximum strain amount of K • ( t /2 ) are generated 
in the opposite directions, wherein the curvature is set to/c 
( = 1/rd; rd being the radius of die shoulder) (Fig. 9(b)). 
Further, here, when an amount of work hardening is ignored (Fig. 
6 explained above) , the distribution of the stress cr is made such 
that the stresses having the maximum stress of CTp are generated 
on the front surface side and the rear surface side in the 
different directions (compression stress and the tensile stress ) 
(Fig. 9(c)). Due to such a stress distribution, the bending 
moment is generated. Then, assuming the position in the sheet 
thickness direction as v and the distance from the sheet 
thickness center to both end portions of the elastic region as 
+yi , -yi ( Fig . 9(c)), yi becomes yi = cTp/ ( E • fc ) wherein E is the 
Young's modulus and /c is the curvature. Accordingly, the 
bending moment acting on the central portion of the metal sheet 
is expressed by a following equation (IV). 



Further, in the stage ( unbending stage ) shown in Fig. 7(c), 
in the marked position X of the material, the whole region once 
returns to the elastic region and thereafter the region is yielded 



2 



2 




(IV) 
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again in the direction opposite to the direction of the stage 
of Fig. 7(b) (Fig. 10(a)) so that the distribution of the strain 
a is not generated (Fig. 10(b) ) . Further, here, when an amount 
of work hardening is ignored, the distribution of the stress G 
is made such that the stresses having the maximum stress of O 
p ( crp being the yield stress) are generated in the different 
directions from each other and in the directions opposite to the 
directions of Fig. 9(c) (Fig. 10(c)). Here, assuming the 
position in the sheet thickness direction as V and the distance 
from the sheet thickness center to both end portions of the 
elastic region ( elastic-plastic boundary) as +y2, -y2 (Fig. 
10(c)), y2 becomes y2 = 2crp/(E-A:), wherein E is the Young's 
modulus and IC is the curvature. Accordingly, the bending moment 
amount M^acting on the central portion of the metal sheet is 
expressed by a following equation (V) . 



Then, at the time of removing the metal sheet from the 



mold, the external force is released and hence, the whole region 
in the sheet thickness direction of the metal sheet becomes the 
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elastic region (Fig. 11(a)) and the elastic recovery (wall warp 
p ) is generated such that the bending moment amount which the 
region has due to the elasticity of the material becomes zero. 
Here, the distributions of the strain and stress which are 
generated in the sheet thickness direction are formed as shown 
in Fig. 11(b) (maximum strain £ = P * (t/2)) and Fig. 11(c). 

Then^ the bending moment amount e reduced by the 
above-mentioned elastic recovery can be expressed by a following 
equation (VI) and this elastic recovery is generated such that 
it cancels the bending moment amount M(3) which is generated by 
the distribution of the stress at the time of forming. 
Accordingly, the wall warp amount p is generated corresponding 
to the bending moment amount M(4) which satisfies a following 
equation (VII ) . 




(VI) 



(VII) 



To put the above-mentioned equation (VII) in order with 
respect to the wall warp amount p , a following equation (VIII) 
is formulated. By putting the relationship of the above- 
mentioned equation (V) into this equation (VIII) and put them 
in order, the above-mentioned equation (I) is formulated. 

p = (12/Et^) • M® (VIII) 
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By the way^r in a usual press forming, the value of 7/3 • 
[(2a'p/Et) • rd]^ in the second term of the right-hand member 
in the above-mentioned equation (I) becomes a value extremely 
smaller than 1. For example, assuming that cTp = 600 MPa, rd = 
5mm, E = 205800 MPa, t = 1.2 mm, the value of 7/3 • [(2crp/Et) • 
rd]^ becomes 1.38 X 10"^ ( < 1 ). Accordingly, the second term 
{ 1 - 7/3 • [(2crp/Et) * rd]^ > of the right-hand member in the 
above-mentioned equation (I) becomes approximately 1 so that it 
can be ignored whereby the equation (I) can be further simplified 
like the above-mentioned equation (II). 

The inventors of the present invention have compared the 
actually measured values of wall warp amount obtained by an 
experiment under various material strengths and forming 
conditions and the predicted values obtained by the above- 
mentioned equation { II ) . The actually measured values used here 
are sampled values at a relatively small region (the region 
surrounded by the broken line in Fig. 4(d) , (e) ) where the radius 
of die shoulder is extremely small and the tension (blank holding 
force) is relatively small. As a result, it was found that, to 
accurately predict the wall warp amount, it is necessary to take 
the work hardening into consideration to some extent. 

Then, to reflect the influence of the work hardening as 
much as possible, the inventors of the present invention have 
made further extensive studies. As a result, it was found that 

-20- 



by setting a value which is equal to or less than the tensile 
stress and exceeds the yield stress CTp as an apparent yield stress 
CTp' which includes an amendment corresponding to the work 
hardening and is different from the reality in place of the yield 
stress CTp of the above-mentioned equation (II), the wall warp 
amount p may be predicted by a following equation (II)' based 
on this apparent yield stress cTp' as shown in Fig. 12. Further, 
it was found that when the above-mentioned apparent yield stress 
CTp' is set to the internal division value obtained by the 
following equation (1) based on the actually measured yield 
strength YS and tensile strength TS or when a coefficient in the 
following equation (1) is set to a function (k = f (TS/t)) of 
the ratio (TS/t) between the tensile strength TS and the sheet 
thickness t, the actually measured value and the predicted value 
relatively favorably agree with each other. 

Here, as the above-mentioned yield strength YS and the 
tensile strength TS, values obtained by a usual tension test may 
be used. Further, the reason that the coefficient k is set to 
the function of TS/t can be explained as follows. In general, 
when the stress-strain diagram is measured by the tension test 
or the like, a tendency of work hardening shown in Fig. 21 is 
observed. With respect to materials having different tensile 
strength and sheet thickness, the degree of work hardening at 
the same strain value (working with same tool) can be summarized 
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as follows. 

(a) Corresponding to the increase of the tensile strength, the 
strain value which reaches TS becomes smaller and the degree of 
work hardening (being evaluated by Y1/Y2 in Fig. 21(a), for 
example) at the same strain value is increased. As a result, 
corresponding to increase of the tensile strength, the yield 
strength after work hardening becomes a value closer to the 
tensile strength TS . 

(b) Corresponding to the decrease of the sheet thickness, the 
strain value which reaches TS becomes smaller and the degree of 
work hardening at the same strain value is increased. As a 
result, corresponding to the increase of the tensile strength, 
the yield strength after work hardening becomes a value closer 
to the tensile strength TS. 

Various forms of functional equations which relate k, TS and 
t with each other in view of these tendencies are considered. 
Among these forms, it was experimentally confirmed that when the 
form which uses k as the function of TS/t is adopted, the actually 
measured values and the predicted values agree with each other. 

p = (3crp'/Et) (II) ' 

where p = wall warp amount (curvature; l/mm),crp': apparent 
yield strength (MPa), E: Young's modulus (MPa) , t: sheet 
thickness (mm) . 

CTp' = k • YS + ( 1 - k) TS ( 1 ) 
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where k is the coefficient of the internal division value. 

Here, in the above-mentioned equation (II) ' , the influence 
of the radius rd of die shoulder is not sufficiently expressed. 
However, when the inventors of the present invention have 
reviewed with respect to the manner of change of the actually 
measured value by an experiment, it was found that when a 
correction is made such that a correction term [-C • (rd - 5)] 
obtained by the experiment is added to the above-mentioned 
equation (II)' to formulate a following equation (II)", the 
prediction accuracy is further enhanced. 

p = (3crp'/Et) -c • (rd - 5) (II)" 
However, even with such an equation (II)", the influence 
of tension (blank holding force) is not sufficiently expressed. 
That is, the equation (II)" is only applicable to a region where 
the tension is relatively small (a region in the vicinity of the 
top of Fig. 4(e)). Accordingly, when the inventors of the 
present invention have reviewed the condition of change of 
actually measured values by an experiment with respect to the 
correction when the tension is large, it was found that when the 
correction which is formulated by a following equation (2) is 
made to the above-mentioned equation (II)", a favorable 
prediction accuracy can be obtained. That is, while satisfying 
the following equation (2) formulated in the above-mentioned 
manner, by setting the internal division value obtained by a 
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following equation (3) (same as the previously mentioned 
equation (1)) based on the actually measured yield strength YS 
and the tensile strength TS obtained by the tension test as the 
apparent yield stress CTp' and by adopting a following equation 
(4) using the coefficient k of the equation (3) as the function 
(k = f (TS/t) ) of the ratio (TS/t) between the tension strength 
TS and the sheet thickness t , the wall warp amount p can be easily 
and accurately predicted. 

p= (3crp'/Et) • { 1 - D • [ ( o-t/TS) - 0.3]^} 

- C • ( rd - 5 ) ( 2 ) 

crp' = k • YS + ( 1 - k)TS (3) 
k = A • (TS/t) + B (4) 
where p = wall warp amount (curvature; l/mm),crp': apparent 
yield strength (MPa), E: Young's modulus (MPa), t: sheet 
thickness (mm) , TS : tensile strength (actually measured value; 
MPa), YS: yield strength (actually measured value; MPa),crT: 
tension acting on wall portion (MPa) , rd: radius of die shoulder 
of press forming tool (mm) , k: coefficient of internal division 
value. A: negative constant, B, C, D: positive constant. 

Subsequently, a case in which the angular change amount 
A 0 is predicted is explained. In the bending forming, a metal 
sheet to be subjected to such a forming receives the bending 
deformation from a shoulder portion or the like of a punch which 
constitutes a press forming tool. Accordingly, at the time of 
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bending forming, in the metal sheet, the difference between 
stresses having opposite signs (the difference between a tensile 
strength and a compression stress) is generated on the front and 
back in the sheet thickness direction in the same manner as the 
deep drawing. It has been known that even when such stress 
difference is balanced at a point where an external force is 
removed at the time of the removal of the metal sheet from the 
mold and hence, the bending moment becomes zero, a part of the 
stress difference remains and generates the above-mentioned 
angular change phenomenon. Further, it has been known that 
following tendencies (1) to (4) exist with respect to the angular 
change amount A ^ . 

(1) Corresponding to the increase of the strength (tensile 
strength TS) of the metal material, the angular change amount 
AO is also increased (Fig. 13(a)). 

(2) Corresponding to the decrease of the sheet thickness 
of the metal sheet, the angular change amount A 0 Is increased 
(Fig. 13(b)). 

(3) Corresponding to the decrease of the Young's modulus 
E of the metal material, the angular change amount A 0 is 
increased (Fig. 13(c)). 

(4) Corresponding to the increase of the radius rd of the 
shoulder of the bending tool (punch or the like), the angular 
change amount A 0 is increased (Fig. 13(d)). 
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To simplify the prediction equation, the inventors of the 
present invention have assumed the stress-strain relationship 
that the material does not generate the work hardening although 
it is different from the actual material behavior, as in the case 
of the above-mentioned prediction equation of the wall warp 
amount p (Fig. 6 mentioned previously) and have tried to 
formulate the simpler prediction equation under such a 
relationship. However, since the angular change defect 
generated at the time of bending forming is generated at portions 
which receive only the bending caused by the shoulder portion 
or the like of the punch, it is unnecessary to take the unbending 
phenomenon into account. 

As a result, it was found that the angular change amount 
A Q after the bending deformation and the elastic recovery can 
be expressed by following simple equations such as an equation 
(5) and an equation (IX). 

1^6= -Q ' (rp + t/2) • A p (5) 
A/0= (-3crp/Et) • [1 - 1/3 (2crp'/Et) -rp]^] (IX) 
where A 0 : angular change amount (degree) , 0 : bending angle 
(degree), rp: radius of shoulder of bending tool (mm), t: sheet 
thickness (mm), A pi curvature change amount (1/mm), cJp': 
apparent yield strength (MPa), E: Young's modulus (MPa). 

The above-mentioned equation (5) and the equation (IX) 
can be formulated in the following manner. That is, the 
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curvature change amount A p of the bending portion due to the 
elastic recovery after bending forming is generated such that 
it cancels the bending moment M(3) as in the case in which the 
previously mentioned equation (VII) is formulated and hence, the 
curvature change amount A p is generated to satisfy a following 
equation (X) . 

M(^-r^E-Ap-ridtj = M^- Ap = 0 (X) 

J- 2 12 

Accordingly, the curvature change amount A p can be 
expressed by a following equation (XI) and the above-mentioned 
equation (IX) can be formulated by putting this equation (XI) 
in order. 

A p = (-12/Et^) • M(2, (XI) 
On the other hand, to consider the geometric conversion of 
the above-mentioned curvature change amount A p , since the length 
of the central portion of the sheet thickness does not change even 
after the curvature is changed, the relationship of a following 
equation (XIII) is established. By putting this equation (XIII) 
in order with respect to the angular change amount A 0 , the 
previously-mentioned equation (5) can be formulated. 
0/[l/[rp + (t/2)]] 

= {d +A(9)/[l/[rp + t/2] - Ap] (XIII) 
where: A0: angular change amount (degree), 0: bending angle 
(degree), rp: radius of shoulder of bending tool (mm), t: sheet 
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thickness (mm), Ap: curvature change amount (1/mm). 

In case a punch having a punch shoulder radius of 
approximately 3-20 mm which is popularly used in a usual press 
forming is used, the above-mentioned equation (IX) can be further 
simplified to a following equation (XII) as in the case of the 
deep drawing. 

A p = (-3crp/Et) (XII) 
Further, in the present invention, since the influence 
of the work hardening is considered in a form as simple as possible, 
as in the case of the prediction equation of the wall warp amount, 
as the value of the yield strength CTp, the values defined by the 
previously mentioned equation (1) are used as the apparent yield 
strength a p' which is obtained by correcting an amount 
corresponding to the work hardening and is different from the 
actually measured value. Further, in the above-mentioned 
equation (XII), since the influence of the shoulder radius of 
the bending tool (punch or the like) is not sufficiently expressed, 
it was found that by judging the manner of changing of the actually 
measured values and performing the correction as expressed by 
a following equation (6), a favorable prediction accuracy can 
be obtained. 

Ap= (-3crp'/Et) • [1 + exp(-G-rp)] (6) 
That is, along with the equation (5) and the equation (6) 
formulated in the above-mentioned manner, by setting the 
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internal division value formulated by a following equation (7) 
(same as the previously mentioned equation (1), (3)) based on 
the actually measured yield strength YS and the tensile strength 
TS obtained by the tension test as the apparent yield stress O 
p' and by adopting a following equation (8) using the coefficient 
k of the equation (7) as the function (k = f (TS/t)) of the ratio 
(TS/t) between the tension strength TS and the sheet thickness 
■t, the angular change amount A 0 can be easily and accurately 
predicted . 

crp' = k • YS + ( 1 - k)TS (7) 
k = A • (TS/t) + B (8) 
where ap': apparent yield strength (MPa), E: Young's modulus 
(MPa), YS: yield strength (actually measured value; MPa),TS: 
tensile strength ( actually measured value; MPa ) , k: coefficient 
of internal division value. A, B, G: constant. 

The advantageous effects of the present invention will 
be further specifically explained in conjunction with 
embodiments hereinafter. However, none of following embodiments 
limits the present invention and any design matters regarding 
the present invention are included in the technical scope of the 
present invention in view of the spirit of the present invention 
which has been explained heretofore and which will be explained 
hereinafter . 
(Embodiment 1) 
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First of all, the prediction accuracy in the region (see 
Fig. 4(e)) where the tension is relatively small was confirmed. 
Using various kinds of steel sheets having the yield strength 
YP, the tensile strength TS and the sheet thickness t shown in 
following Table 1 and under forming conditions (die shoulder 
diameter rd) shown in the Table 1, a "hat channel" specimen 
forming experiment was carried out. Here, BHF (blank holding 
force) was adjusted such that the influence which the tension 
cTt affects the wall warp amount can be made as small as possible 
and the ratio ( cTt /TS) between the tens ion cTt and the tension 
strength TS substantially becomes constant (0.3 i the region in 
the vicinity of the top point in Fig. 4(e)). 



Table 1 





Yield strength 


Tensile strength 


Sheet thickness 


Die shoulder 


Actually 


Predicted value 


No. 


(MPa) 


(MPa) 


(mm) 




measured value 


Peal 








(mm) 


[x10->m-^)] 


[x10>m-')] 


1 


608 


1009 


1.2 


5 


12.196 


12.255 


2 


373 


608 


1.4 


5 


5.100 


5.365 


3 


406 


596 


1.2 


5 


5.788 


6.469 


4 


362 


457 


1.0 


5 


5.393 


6.146 


5 


351 


450 


1.0 


5 


5.571 


6.012 


6 


277 


444 


1.6 


5 


2.407 


3.213 


7 


177 


306 


1.2 


5 


2.048 


2.826 


8 


406 


596 


1.2 


5 


6.697 


6.469 


9 


463 


611 


1.4 


5 


5.151 


5.757 


10 


324 


616 


1.4 


5 


5.177 


5.231 


11 


373 


608 


1.4 


15 


4.390 


4.284 


12 


406 


596 


1.2 


15 


5.691 


5.388 


13 


362 


457 


1.0 


15 


4.172 


5.065 


14 


351 


450 


1.0 


15 


5.634 


4.931 


15 


277 


444 


1.6 


15 


2.483 


2.132 


16 


177 


306 


1.2 


15 


2.565 


1.745 


17 


463 


611 


1.4 


15 


4.398 


4.676 


18 


324 


616 


1.4 


15 


4.172 


4.150 



with respect to the wall warp amount ( curvature ) p in the 



test, the actually measured values { pmes) were shown in the 
Table 1. Using these actually measured values, the values of 
the constants A, B and C in the previously mentioned equations 
(2) to (4) were determined. Here, since the ratio cTt /TS is set 
toc7T/TS = 0.3, the value of D • [ ( cTt/TS) - 0.3]^ does not depend 
on the value of D and becomes zero. Accordingly, the 
consideration is made ignoring D. The method for determining 
the value of D will be explained later. Under such conditions, 
the values of A, B and C were determined using Newton-Raphson ' s 
method which is one of optimization methods such that the square 
values of the differences between the actually measured values 
and the predicted values by these equations can be minimized. 
As a result, the values of A, B and C were respectively set to 
A = -9.7 08 X 10"^ (mm/MPa) , B = 0.8161,C = 1.082 X 10"^ ( 1/mm^) . 
The obtained values of A,B and C and the predicted values (p 
cal) obtained using the equations (2) to (4) are also shown in 
the Table 1. The result obtained by comparing the actually 
measured values (pmes) and the predicted values (peal) is shown 
in Fig. 14. The high correlation was recognized between the 
actually measured values (pmes) and the predicted values (p 
cal) so that it is understood that the wall warp amount can be 
predicted with a favorable accuracy. 

Then, it was confirmed whether the above-mentioned 
equations (2) to (4) and the values of A, B and C obtained by 
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the above-mentioned experiment are applicable to the forming of 
aluminum sheets which largely differ in the Young's modulus from 
the steel sheets used in the above-mentioned experiment or not. 
Using various kinds of aluminum sheets (3000 series, 7000 series) 
having the yield strength YP, the tensile strength TS and the 
sheet thickness t shown in following Table 2 and under forming 
conditions (die shoulder diameter rd) shown in the Table 2 and 
setting the ratio (^t /TS tocTT/TS = 0.3 (see Fig. 4(e) previously 
mentioned), a "hat channel" specimen forming experiment was 
carried out. In the same manner as the experiment on the steel 
sheets, the actually measured values (pmes) and the predicted 
values (pcal) based on the previously mentioned equation (2) 
to (4) of the wall warp amount p (curvature) were compared with 
each other and the result is shown in Fig. 15. The high 
correlation was recognized between the actually measured values 
( pmes) and the predicted values ( /Ocal) so that it is understood 
that the wall warp amount can be predicted with a favorable 
accuracy. 



Table 2 





Yield strength 


Tensile strength 


Sheet thickness 


Die shoulder 


Actually measured 


Predicted value 


No. 


(MPa) 


(MPa) 




(mm) 


value 
Pn,es[x10">m-')] 


Pcal 

[XlO->m-^)l 


19 








3 


5.184 


5.018 


20 


73 


177 


1.0 


5 


5.039,4.777 


4.802 


21 








10 


4.202 


4.261 


22 


58 


153 


1.0 


5 


5.049 


3.910 


23 


131 


284 


1.0 


5 


9.965 


8.817 



Subsequently, the prediction accuracy in the region in 
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which the tension (Tt is relatively large was also confirmed . Here, 
the adjustment of the coefficient D was conducted so as to correct 
the influence of the tension. Here, the actually measured values 
of the wall warp amount were obtained with respect to the specimen 
No. 2 in the Table 1 by changing the tension CTt to 200 MPa, 300 
MPa and 450 MPa with the adjustment of the BHF (blank holding 
force) . Using these actually measured values and the values of 
the constants A, B and C obtained by the previous experiment, 
the value of D was determined using Newton-Raphson ' s method such 
that the square values of the differences between the actually 
measured values and the predicted values obtained by the 
equations (2) to (4) can be minimized. As a result, D was set 
to D = 3.559. The result obtained by comparing these actually 
measured values ( pmes) and the predicted values ( peal) is shown 
in Fig. 16. The high correlation was recognized between the 
actually measured values (pmes) and the predicted values (p 
cal) so that it is understood that the wall warp amount can be 
predicted with a favorable accuracy. 
(Embodiment 2) 

With respect to the prediction accuracy of the prediction 
equations of the angular change amount A 0 , the evaluation was 
performed using the basic bending test. First of all, using 
various kinds of steel sheets having the yield strength YP, the 
tensile strength TS and the sheet thickness t shown in following 
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Table 3 and under forming conditions (punch shoulder radius rp) 
shown in the Table 3, an L-bending forming experiment was carried 
out. The state of the L-bending forming experiment is shown in 
Fig. 17. With respect to the angular change amount A 0 in the 
test, the actually measured values (A^mes) were shown in the 
Table 3. 
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Table 3 



No. 


Yield strength 
(MPa) 


Tensile strength 
(MPa) 


Sheet thickness 
(mm) 


Die shoulder 
radius (mm) 


Actually 
measured value 

Pmes 

[x10->m-')] 


Predicted value 

Peal 

[XlO-^(mm-')] 


24 
25 
26 


176 


318 


1.2 


1 

5 
10 


0.95 
2.10 
2.70 


0.795 
2.342 
3.771 


27 
28 
29 


425 


647 


1.2 


1 
5 
10 


4.25 
5.00 
8.75 


1.938 
5.712 
9.197 


30 
31 
32 


672 


1032 


1.2 


1 
5 

10 


4.70 
10.30 
18.15 


3.437 
10.128 
16.307 


33 
34 
35 


1183 


1327 


1.2 


1 
5 
10 


5.65 
12.90 
19.65 


4.419 
13.023 
20.969 


36 
37 
38 


373 


608 


1.4 


1 
5 
10 


2.55 
5.00 
7.65 


1.563 
4.412 
7.045 



Using these actually measured values, the constants A, 



B and C in the previously mentioned equations (5) to (8) were 
determined. In the same manner as the prediction equation of 
the wall warp amount p , the constants A, B and C are coefficients 
provided for determining the apparent yield stress (7 p'. 
Accordingly, these values were set identical with the values 
obtained in the previously mentioned experiment for the wall warp 
amount. Here, only the coefficient G which is provided for 
correcting the influence of the punch shoulder radius rp was 
determined using Newton-Raphson ' s method such that the square 
values of the differences between the actually measured values 
and the predicted values obtained by these equations can be 
minimized. As a result, the value of G was set to G = 0.1012 
(mm) . The obtained coefficient values and the predicted values 
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obtained using the equations (5) to (8) are also shown in the 
Table 3. The result obtained by comparing the actually measured 
values (A^mes) and the predicted values (A^cal) is shown in 
Fig. 18. The high correlation was recognized between the 
actually measured values (A^mes) and the predicted values (A 
^cal) so that it is understood that the angular change amount 
can be predicted with a favorable accuracy. 

Subsequently, to confirm the prediction accuracy in a 
region where the punch shoulder radius rp is relatively large, 
using steel sheets having the yield strength YP, the tensile 
strength TS and the sheet thickness t shown in following Table 
4 and under forming conditions (punch shoulder radius rp) shown 
in the Table 4, a U-bending forming experiment was carried out. 
The state of the U-bending forming experiment is shown in Fig. 
19 . 



Table 4 





Yield strength 


Tensile strength 


Sheet thickness 


Die shoulder 


Actually 


Predicted value 


No. 


(MPa) 


(MPa) 


(mm) 


radius (mm) 


measured value 

Pnes 

[xlO-^mm-')] 


[x10-=(mm-')] 


39 








4.8 


8.987 


8.619 


40 


375 


610 


1.4 


9.0 


15.063 


13.195 


41 








18.0 


20.886 


21.076 


42 








38.0 


38.734 


38.346 



The prediction of the angular change amount A d was 
carried out using the previously mentioned values of A, B, C and 
G and the previously mentioned equations (5) to (8) . The result 
obtained by comparing the actually measured values (A 0mes) and 
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the predicted values (A^cal) obtained from the previously 
mentioned equations (5) to (8) is shown in Fig. 20. The high 
correlation was recognized between the actually measured values 
(A^mes) and the predicted values (A^cal) so that it is 
understood that the angular change amount can be predicted with 
a favorable accuracy. 

According to the present invention, with the 
constitutions which have been explained heretofore, even when 
engineers having no experience and technical storage do not have 
expertise such as numerical simulations, they can easily and 
accurately predict an amount of dimensional accuracy defect at 
the time of press-forming. Such a method enables engineers to 
take a prompt and effective measures to enhance the dimensional 
accuracy against the dimensional accuracy defect which has been 
considered as a particularly important problem recently and 
hence the expectation on the technical significance of this 
method is considered enormous. 
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